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Abstract. Let Slbea bounded domain with C 2 -smooth boundary in an n-dimensional 
oriented Riemannian manifold. It is well-known that for the bi-harmonic equation 
A 2 m = in f! with the O-Dirichlet boundary condition, there exists an infinite set {u^} of 
biharmonic functions in Q with positive eigenvalues {A^} satisfying Ati^ + XkB^^r = 
on the boundary c9f2. In this paper, by a new method we establish the Weyl-type 
asymptotic formula for the counting function of the biharmonic Steklov eigenvalues . 



1. Introduction 

Spectral asymptotics for partial differential operators have been the subject of extensive 
research for over a century. It has attracted the attention of many outstanding mathemati- 
cians and physicists. Beyond the beautiful asymptotic formulas that are intimately related 
to the geometric properties of the domain and its boundary, a sustaining force has been 
its important role in mathematics, mechanics and theoretical physics (see, for example, [4], 
0, i, 0, i, [18], m US, E3, [32], [33], [36], [39], [40], [41], 02], [50]). 

Let (M,g) be an oriented Riemannian manifold of dimension n with a positive definite 
metric tensor g, and let O C M. be a bounded domain with C 2 -smooth boundary 90. 
Assume g is a non- negative bounded function defined on 90. We consider the following 
classical biharmonic Steklov eigenvalue problem: 

A 2 g u = in O, 

(1.1) { u = on 90, 



A g u + Xgp = on 90, 



where v denotes the inward unit normal vector to 90, and A g is the Laplace-Beltrami 
operator defined in local coordinates by the expression, 
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Here \g\ := det(gij) is the determinant of the metric tensor, and g n are the components of 
the inverse of the metric tensor g. 

The problem (jl.ip has nontrivial solutions u only for a discrete set of A = A& , which are 
called biharmonic Steklov eigenvalues (see [TT], [2JJ, [35] or [17]). Let us enumerate the 
eigenvalues in increasing order: 

< Ai < A 2 < • • • < A fc < • • • , 

where each eigenvalue is counted as many times as its multiplicity. The corresponding 
eigenfunctions ^7, ■ • • , ^7, • • ■ form a complete orthonormal basis in L 2 JdVL) (see, 
Proposition 3.5). It is clear that Xk can be characterized variationally as 

J n \A gUl \ 2 dR . /„ \A g v\ 2 dR 

Al = 5 = lilt 5 , 

Afe = = — = max mi 5 — , fc = 2,3,4, ••• 

where H m (il) is the Sobolev space, and where dR and ds are the Riemannian elements of 
volume and area on and <9f2, respectively. 

In elastic mechanics, when the weight of the body ft is the only body force, the stress 
function u must satisfy the equation A 2 it = in 51 (see, p. 32 of [45 ). In addition, the 
boundary condition in (jl.ll) has an interesting interpretation in theory of elasticity. Con- 
sider the model problem (see [TT]): 

A 2 u = f in tt, 

u = 0, Ait + (1 - a)i = on dfl, 

where 57 C M 2 is an open bounded domain with smooth boundary, a G ( — 1, 1/2) is the 
Poisson ratio and 1 is the mean curvature of the boundary dfl. Problem (|1.2|) describes the 
deformation u of the linear elastic supported plate fl under the action of the transversal 
exterior force / = /(x), x G Q. The Poisson ratio a of an elastic material is the negative 
transverse strain divided by the axial strain in the direction of the stretching force. In 
other words, this parameter measures the transverse expansion (respectively, contraction) 
if er > (respectively, a < 0) when the material is compressed by an external force. We refer 
to [22], [46] for more details. The restriction on the Poisson ratio is due to thermodynamic 
considerations of strain energy in the theory of elasticity. As shown in [22], there exist 
materials for which the Poisson ratio is negative and the limit case a = — 1 corresponds 
to materials with an infinite flexural rigidity (see, p. 456 of [41]). This limit value for a is 
strictly related to the eigenvalue problem Hence, the limit value a = — 1, which is not 

allowed from a physical point of view, also changes the structure of the stationary problem 
(|1.2p : For example (see [JT]), when fi is the unit disk and Ai = (1 — 0)1 = 1 — a = 2, (jl.2| 
either admits an infinite number of solutions or it admits no solutions at all, depending on 
/■ 

Problem (jl.ip is also important in conductivity and biharmonic analysis because the 
related problem was initially studied by Calderon (cf. [3j). This connection arises because 
the set of the eigenvalues for the biharmonic Steklov problem is the same as the set of 
eigenvalues of the well-known "Neumann-to-Laplacian" map for biharmonic equation (This 
map associates each normal derivative du/dv defined on the boundary dtt to the restriction 
(Au)| afJ of the Laplacian of u for the biharmonic function u on O, where the biharmonic 
function u is uniquely determined by u\ gQ = and (du/dv) | ao )- 



(1.2) 
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In the general case the eigenvalues A& can not be evaluated explicitly. In particular, for 
large k it is difficult to calculate them numerically. In view of the important applications, 
one is interested in finding the asymptotic formulas for as k — > oo. However, for a 
number of reasons it is traditional in such problems to deal with the matter the other way 
round, i.e., to study the sequential number k as a function of r. Namely, let us introduce 
the counting function A(r) defined as the number of eigenvalues less than or equal to 
a given r. Then our asymptotic problem is reformulated as the study of the asymptotic 
behavior of A(t) as t — > +00. 

In order to better understand our problem and its asymptotic behavior, let us 

mention the Steklov eigenvalue problem for the harmonic equation 



A g v = in 17, 

dv ' rjgv = on dfl, 



( L3 ) 

dv 

where 77 is a real number. This problem was first introduced by V. A. Steklov for bounded 
domains in the plane in [41] (The reader should be aware that "Steklov" is also often 
transliterated as "Stekloff".) His motivation came from physics. The function v represents 
the steady state temperature on f2 such that the flux on the boundary is proportional to the 
temperature (In two dimensions, it can also be interpreted as a membrane with whole mass 
concentrated on the boundary). For the harmonic Steklov eigenvalue problem (|1.3[) . in a 
special case in two dimensions, A. Pleijel [35] outlined an investigation of the asymptotic 
behavior of both eigenvalues and the eigenfunctions. In 1955, L. Sandgren [35] established 
the asymptotic formula of the counting function B(t) = #{r//c|% < t}: 

(1.4) B(t) ~ "T\ -. / g^ds asT^+oo, 



(2?r) n - 



i.e. 



hm — -^r = -— — - / g ds, 

where w n _i is the volume of the unit ball of R ra_1 , and the integral is over the boundary 
dtt. This asymptotic behavior is motivated by the similar one for the eigenvalues of the 
Dirichlet Laplacian. The classical result for the Dirichlet (or Neumann) eigenvalues of the 
Laplacian on a smooth bounded domain is Weyl's formula (see [4], [6] or [49]): 

(1.5) N(t,Q) ~ ^(vol(ft))W 2 as r ^ +00, 

where N(r, fl) is the number of the Dirichlet (or Neumann) eigenvalues of domain f2 less 
than or equal to a given r. In the case of two-dimensional Euclidean space, Pleijel [34] in 
1950 proved an asymptotic formula for the eigenvalues ~E? k of a clamped plate problem: 

A 2 u - E 2 u = in Q, 

u = ^ = o on an. 

Of 

Grub [IB] and Ashbaugh, Gesztesy, Mitrea and Teschl [2] obtained Weyl's asymptotic for- 
mula for the eigenvalues A& of the buckling problem in R n : 

\ 2 u + AAu = in 0, 
, = |2i = on 9S1. 



(1.6) 



(1.7) 



Note that for the Dirichlet eigenvalues, the Neumann eigenvalues, the buckling eigenvalues 
and the square root of the clamped plate eigenvalues in a fixed domain, their counting 
functions have the same asymptotic formula (|1.5p (see, for example, @], [3S], [13], [35] and 



4 



GENQIAN LIU 



The study of asymptotic behavior of the biharmonic Steklov eigenvalues is much more 
difficult than that of the harmonic Steklov eigenvalues. It had been a challenging problem 
in the past 50 years. The main stumbling block that lies in the way is the estimates for the 
distribution of the boundary eigenvalues for bi-harmonic equations with suitable boundary 
conditions. Some important works have contributed to the research of this problem, for 
example, L. E. Payne 32 , J. R. Kuttler and V. G. Sigillito [21], A. Ferrero, F. Gazzola 
and T. Weth [H, Q. Wang and C. Xia }47], and others. 

In this paper, by a new method we establish the Weyl-type asymptotic formula for the 
counting function of the biharmonic Steklov eigenvalues. The main result is the following: 

Theorem 1.1. Let (A4,g) be an n- dimensional oriented Riemannian manifold, and let 
£1 C M. be a bounded domain with C 2 -smooth boundary dfl. Then 



where A(t) is defined as before. 

Corollary 1.2. Let (M.,g) be an n-dimensional oriented Riemannian manifold, and let 
f2 C M. be a bounded domain with C 2 -smooth boundary dfl. If, in problem U.l\) , g = 1 on 
dQ, then 



We outline the idea of the proof of Theorem 1.1. First, we make a division of f2 into 
subdomains (by dividing dQ into sufficiently small parts, then taking a depth a > (small 
enough) in the direction of inner normal of dQ to form a finite number of n-dimensional 
subdomains). Under a sufficiently fine division of d£l (also a sufficiently small), g lk and g 
can be replaced by constants because their variant will be small, so that the corresponding 
subdomains whose partial boundaries are situated at the 90 can be approximated by 
Euclidean cylinders. Next, we construct three Hilbert spaces of functions and their self- 
adjoint linear transformations whose eigenfunctions are just the Steklov eigenfunctions 
with corresponding boundary conditions. It can be shown that these Steklov eigenvalue 
problems have the same boundary conditions on the base of each cylinder as the original 
one in problem (jl.lj) but they have relevant boundary conditions on the other parts of a 
cylinder. In particular, on each cylindrical surface, these boundary conditions will be one of 
the three forms u = A g u — 0, A g u = d ^l^ = and |^ = — 0. The main purpose 

of constructing such Steklov problems is that when putting together such cylinders, we 
can obtain global upper and lower estimates for the counting function A(t) of the original 
Steklov problem (i.e., A°(t) < A(t) < A d (r) < A f (t) for all r > 0, see Sections 3, 6). For 
each Euclidean cylinder, by using a cubical net we can divide the base of the cylinder into 
(n — l)-dimensional cubes and some smaller parts which intersect boundary of the base, so 
that we get n-dimensional parallelepipeds and some smaller n-dimensional cylinders. As for 
the n-dimensional parallelepiped, we can explicitly calculate the Steklov eigenfunctions and 
eigenvalues by separating variables, and then we can compute the asymptotic distribution 
of eigenvalues by means of the well-known variational methods used by H. Weyl [50], R. 
Courant and D. Hilbert [6 in the case of the membranes. Meanwhile, for each small n- 
dimensional cylinder, by introducing a nice transformation we may map it into a special 
cylinder whose counting functions of Steklov eigenvalues can also be estimated. Finally, 
applying normal coordinates system at a fixed point of each subdomain of a division and 



(1.8) 





(1.9) 
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combining these estimates, we establish the desired asymptotic formula for A(t). Note that 
the Holmgren uniqueness theorem for the solutions of elliptic equations plays a crucial role 
in this paper. 

This paper is organized as follows. In Section 2, we prove two compact trace lemmas 
for bounded domains with piecewise smooth boundaries. In Section 3, we define various 
self-adjoint transformations on the associated Hilbert spaces of functions, and give the 
connections between the eigenfunctions of self-adjoint transformations and the Steklov 
eigenfunctions (corresponding to different kinds of boundary conditions). Section 4 is 
dedicated to deriving the explicit formulas for the biharmonic Steklov eigenvalues and 
eigenfunctions in an n-dimensional rectangular parallelepiped of R™ , which depends on a 
key calculation for the solutions of biharmonic equations. The counting functions of Steklov 
eigenvalues for general cylinder of the Euclidean space are dealt with in Section 5. In the 
final section, we prove Theorem 1.1 and Corollary 1.2 on Ricmannian manifolds. 



An n-dimensional cube in R n is the set {x £ R™|0 < Xi < a, i = 1, • ■ • , n}. 

Let / be a real- valued function defined in an open set fi in R™ (n > 1). For y G fi we 
call / real analytic at y if there exist ap G M 1 and a neighborhood U of y (all depending 
on y) such that 



for all x in U. We say / is real analytic in f2, if / is real analytic at each 

Let O together with its boundary be transformed pointwise into the domain fl' together 
with its boundary by equations of the form 



where the functions fa and their first order derivatives are Lipschitz continuous throughout 
the domain, and they are less in absolute value than a small positive number e. Then we 
say that the domain fl is approximated by the domain fi' with the degree of accuracy e. 

Let (M,g) be a Riemannian manifold. A subset T of (M,g) is said to be an (n — 1)- 
dimcnsional smooth (respectively, real analytic) surface if T is nonempty and if for every 
point x in T, there is a smooth (respectively, real analytic) diffcomorphism of the open unit 
ball 5(0, 1) in R" onto an open neighborhood U of x such that 5(0, l)n{i£ R n \x n = 0} 
maps onto U R T. 

An (n — l)-dimensional surface T in (M,g) is said to be piecewise smooth (respectively, 
piecewise real analytic) if there exist a finite number of (n — 2)-dimensional smooth surfaces, 
by which T can be divided into a finite number of (n— l)-dimensional smooth (respectively, 
real analytic) surfaces. 

A subset of L 2 (T) is called precompact if any infinite sequence of elements of £ 
contains a Cauchy subsequence {uk>} 7 i.e., one for which 



2. Compact trace Lemmas 



/O) = ^a {x-y) 



(2.1) 



x n ), i — 1,2, ••• 



(2.2) 




as fc', I' — > oo. 
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From here up to Section 5, let M. be an n-dimensional Riemannian manifold with real 
analytic metric tensor g. 

Lemma 2.1. Let D C (A4,g) be a bounded domain with piecewise smooth boundary. 
Assume that 9Jt is a set of functions u in Hq(D) D H 2 (D) for which 

(2.3) / \A g u\ 2 dR 

Jd 

is uniformly bounded. Then the set {f^ju 6 931} is precompact in L 2 {dD). 
Proof. Put 

, s , , fn \& Q u\ 2 dR 

2.4 Ai(D) = inf J P 9 ' . 

ueHi(D)nH 2 (D) J D \\7gu\ 2 dR 

We claim that Ai(_D) > 0. In fact, by applying Green's formula (see, for example, [4] or 
[39] ) and Schwarz's inequality we see that for any u e Hq(D) Pi H 2 (D), 



\V g u\ 2 dR 



-u(A g u)dR 



( I u 2 dRj \A g u\ 2 dR 



i.e., 



where 



S D \u\UR - S D \W,u\W 



Since the first Dirichlet eigenvalue Ai(-D) is positive for the bounded domain D, i.e. 

f n \V„u\ 2 dR 

(2.6) 0<A 1 ( J D)= inf Joi 9 J , 

ueiJoMo) J D \u\ 2 dR 

we find by (|2.5|) and (|2.6p that Ai(D) > 0, and the claim is proved. 



From (j276l) and (12741) we obtain that 
(2.7) 



/ \u\ 2 dR<—^— [ \\7 g u\ 2 dR for all u G H^(D) 
Jd Jd 



and 

(2.8) / \\7 g u\ 2 dR< — [ \A g u\ 2 dR for all u E Hq(-D) n H 2 (D). 

Since 9Z? is piecewise smooth, we can write 3D = LI^Ti, where Ti is an (n — 1)- 
dimensional surface. For each fixed i, (i = 1, • • • ,m), we choose a smooth (n — 1)- 
dimensional surface f CC f such that dT' { — dTi and Ti U bounds an n-dimensional 
Lipschitz domain D[ satisfying D[ CC D\JT % . Note that w = on T % for u e H^(D)nH 2 (D) 
(see, for example, p. 62 of [24] or Corollary 6.2.43 of [16]). It follows from the a priori es- 
timate of the elliptic operators (see, for example, Theorem 9.13 of [13]) that there exists a 
constant C; > depending only on n, Ti, D[ and Z? such that 

(2-9) ll«l|fl»(D'0 < a(||Au|| L 2 (D) + |M| L 2 (I)) ). 

By assumption, we have J D \Au\ 2 dR < C for all u G 9Jt, where C > is a constant. 
According to (1277) . (j2~8l) and (|2~9l) . we see that for every »6f, 

(2-10) ||«||fli( D {) < CJ', 
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where Cf > is a constant depending only on n,Fi,D-, D and C. Since D\ is a domain 
with Lipschitz boundary in (M.,g), it follows from the Neumann trace theorem (see, p. 16 
of [2J, p. 127 of [2J], pEU or Chs V, VI of 0) that 



dv 



is precompact for each i (i = 1, • • • , m). Consequently, we obtain that {f^|u € is 
precompact in L 2 (dD). □ 

Lemma 2.2. Lei [M.,g) be a real analytic Riemannian manifold, and let D C (A4,g) 
be a bounded domain with piecewise smooth boundary. Suppose T± is a domain in dD 
with dD — Ti^0 and assume that T 2 is an (n — 1)- dimensional real analytic surface in 
dD satisfying T 2 CC dD — T%. Assume € is a set of functions u in K d (D) = {u\u £ 
H 2 (D), u = on Ti, u = |f = on T 2 } for which 

(2.11) / \A g u\ 2 dR 



<j' 

is uniformly bounded. Then the set {§f | r : u G (£} is precompact in L 2 (Ti). 

Proof. Since dD is piecewise smooth, it follows that Ti can be divided into a finite 
number of smooth (n — 1) dimensional surfaces. Without loss of generality, we let Ti itself 
be a smooth (n — 1) dimensional surface. Put 



(2.12) A ri (£>)= inf 



Id \^A 2 dR 



veK d (D), f D \v\ 2 dR=l J D \v\ 2 dR 

In order to prove the existence of a minimizer to (|2.12p . consider a minimizing sequence 
v m in K d (D), i.e., 

/ \A g v m \ 2 dR -> A Fl {D) =0 as m ->■ +00 

JD 

with J D \v m \ 2 dR = 1. Then, there is a Constance C > such that 

(2.13) \\A g v m \\ L 2 (D ) <C, \\v m \\ L 2 {D ) < C for all m. 

Let {Di} be a sequence of Lipschitz domains such that D\ C D 2 C • • • C -D; C • • • CC 
D UTi UT 2 , Ug-L-D; = D, and Ti UT 2 C dDi for all I. It follows from the a priori estimate 
for elliptic equations (see, for example, Theorem 9.13 of \l'S\ ) that there exists a constant 
C[ > depending only on n, Di, D, Ti and r 2 , such that 

(2.14) \\Vm\\rmDi) < C'l(\\AgV m \\L 2 (D) + \\v m \\L*(D))- 

From this and (|2.13p . we see that 

||Um||i? 2 (A) ^ C i for a11 m -> 

where C" is a constant depending only on n, Di,D, Ti, T 2 and C. For each /, by the Banach- 
Alaoglu theorem we can then extract a subsequence {f;, m }m=i 01 {"m}, which converges 
weakly in H 2 (Di) to a limit u, and converges strongly in L 2 (Di) to u. We may assume 
that {vi+i, m }^ =1 is a subsequence of {vi im }^ =1 for every L Then, the diagonal sequence 
{ v i,i}fZi converges weakly in H 2 to u, and strongly converges to u in L 2 , in every compact 
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subset E of D. It is obvious that ||m||l 2 (_d) = 1- Since the functional J D \A g u\ 2 dR is lower 
semicontinuous in the weak H 2 (Di) topology, we have 

/ \A g u\ 2 dR< lim / \A g v k ,k\ 2 dR, 

J Di fc->oo J Di 

I \A g u\ 2 dR = lim / \A g u\ 2 dR < lim ( lim / | A g v k , k \ 2 dR 
lim (lim / \A g v k , k \ 2 dR) = \rA D )- 



so that 



i- 

T2l 



For each fixed I, since v k ^ k — > u weakly in H (Di), we get that v k ^ k — > u strongly in H r {Di 

dv 



for any < r < 2. Note that = and v k J = for i = 1,2. It follows that 



= and uL = i^L =0. Therefore u E K d (D) is a minimizcr. 

I 1 2 CI' II o v ' 



We claim that X ri (D) > 0. Suppose by contradiction that A ri (D) = = 0. 

Then A 9 u = in D. Since the coefficients of the Laplacian are real analytic in D, and since 
1^2 is a real analytic surface, we find with the aid of the regularity for elliptic equations 
(see, Theorem A of [31], [30] or pQ) that u is real analytic up to the partial boundary IV 
Note that u — ^ ~ on T 2 . Applying Holmgren's uniqueness theorem (see, Corollary 5 
of p. 39 in [37] or p. 433 of [33]) for the real analytic elliptic equation A g u — in D, we get 
u = in D. This contradicts the fact J D \u\ 2 dR — 1, and the claim is proved. Therefore 
we have 

(2.15) / \u\ 2 dR < - 1 / \A g u\ 2 dR for u E K d (D). 
Jd *tA d )Jd 

According to the assumption, there is a constant C" such that 

(2.16) \\A g u\\ L 2 [D) < C" for all u E £. 

Again, applying the a priori estimate for the elliptic equations in some (fixed) subdomain 
Di CCfUTiU T 2 (see, Theorem 9.13 of [13]), we obtain that 

(2-17) IMIiP(A) < C[{\\A g u\\ L 2 (D) + \\u\\ L 2 {D) ), 

where the constant C[ is as in (|2.14[) . By (|2 . 15[) — (|2.17p . we get that for every u E <£, 

Mh HDi) < C", 

where C" > is a constant depending only on n, Di, D, Ti, T 2 and C" . It follows from the 
Neumann trace theorem (see, p. 16 of [2], [14] or [27] ) that {f^| r is precompact 

ini 2 (ri). □ 

The following two results will be needed later: 

Proposition 2.3 (see, p. 12 of [33]). Let II be an isometric transformation which 
maps a Hilbert space H onto a subspace n°7^° of another Hilbert space T-L, so that 

(u°, v°)° = (Il u a , n%°) for all u°, v° E H°. 

Suppose that G° and G are two non-negative, self-adjoint, completely continuous linear 
transformations on H and T-L respectively, such that 

(G°u°, v°)° = (Gn°u°, IlV) for all u°, v° E H°. 
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Then 

A 1 " < Mft for k= 1,2,3,- •• , 
where and {/ifc} are the eigenvalues of G° and G, respectively. 

Proposition 2.4 (see, p. 13 of [39]). Assume thatT-L is a direct sum of p Hilbert spaces 

and that the self-adjoint, completely continuous linear transformation G maps every T-Lj 
into itself, 

GUjdUj, j = 1,2,3,- ■■ ,p. 

Denote by Gj the restriction of G to Hj. Then the set of eigenvalues of the transformation 
G (each eigenvalue repeated according to its multiplicity) is identical to the union of the 
sets of eigenvalues of G±, ■ ■ ■ , D p . 



3. Completely continuous transformations and eigenvalues 

Let (J\4,g) be an n-dimensional real analytic Riemannian manifold and let D C M be 
a bounded domain with piecewise smooth boundary T. Suppose that g is a non- negative 
bounded function defined on T or only on a portion T e of T (measure T g = J r ds > 0) and 
assume that J r gds > 0. In case T e ^ T we denote To = T — T e , and assume that Too is 
a real analytic (n — l)-dimensional surface in IV 

If T e 7^ T (measure T > 0), we denote 

K(D) = {u\u e H%(D) n H 2 (D), and % = on r 00 } 
K d {D) = {u\u e H 2 (D), u = on T e , and u = = on r 00 }. 

If r e = T, we denote 

N(D) = {u\u e Hl{D) n H 2 (D)}. 

It follows from the property of Hq(H) (see, for example, p. 62 of 24 or Corollary 6.2.43 
of [TB] or [29J that u = on 3D for any u E Hq(CI) (Therefore, we always have that u = 
on T for any u E K{D) or N{D)). 

We shall also use the notation 

(u,v)* = [ (Agu)(A g v)dR, u,veK(D) or K d (D) or N(D). 
Jd 

The bilinear functional (u, v}* can be used as an inner product in each of the spaces K(D), 
K d {D) and N(D). In fact, (u, v)* is a positive, symmetric, bilinear functional. In addition, 
if (u, u)* = 0, then A g u = in D. In the case u £ K(D) or N(D), by applying the 
maximum principle, we have u = in D. In the case u G K d (D), since u — ^ = on r 00 , 
we find by Holmgren's uniqueness theorem (see, Corollary 5 of p. 39 in [37]) that it = in 
D. Closing K(D), K d (D) and N{D) with respect to the norm ||it||* = y/(u,u)*, we get 
the Hilbert spaces (fC, \\ ■ ||*), (K d , || ■ ||*) and (TV, || • ||*), respectively. 
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Next, we consider two linear functionals 

du dv 



[u, v] = 



dv dv 



and 

(3.1) (u,v) = (u,v)* + [u,v], 

where u,v E K(D) or u,v G K d (D) or u,v G N(D). It is clear that (u,v) is an inner 
product in each of the spaces K(D), K d (D) and N(D). 

Lemma 3.1. The norm 



\\u\\* = \J (u, u}* 

and 



\\ U \\ = y/(u,u) 

are equivalent in K(D), K d (D) and N(D). 

Proof. Obviously, < ||u|| for all u in K(D) or K d {D) or N(D). In order to prove 
the equivalence of the two norms, we first consider the case in linear space N{D). It suffices 
to show that ||u|| is bounded when u belongs to the set 

m = {u\u G N(D),\\u\\* < 1}. 

It follows from Lemma 2.1 that Wlr ■= {f^| u S Wl} is precompact in L 2 (r). This implies 
that there exists a constant C > such that J r (|^) 2 ds < C for all u G 371. Therefore, 

[u, u] — f r g (§f) 2 ds is bounded in 971, and so is ||w,|| 2 = (u, u)* + [u, u]. Similarly, applying 
Lemmas 2.1, 2.2 we can prove the corresponding results for the spaces K(D) and K d (D). 
□ 

From Lemmas 2.1, 2.2, it follows that 



I [u, u] 



g[ dv ] dS 



<C{u,u)* for all u in K{D) or K d (D) or N{D). 



Therefore, [u,v] is a bounded, symmetric, bilinear functional in (K(D), (•,•)*), (K d (D), (•,•)*) 
and (N(D), (■,■>*). Since it is densely defined in (£,(■,■)*), (IC d , (;■)*) and (A/ - , (■, •}*), 
respectively, it can immediately be extended to (£,(•, •>*), (/C d ,(-, •)*) and (TV, (•, •>*). We 
still use [u, v] to express the extended functional. Then there is a bounded linear transfor- 
mation G ( * } of (AC, ■)*) into (/C, (•, ■)*) (respectively, G ( *l of (/C d , (•, ■)*) into (JC d , (•, ■>*), 
of (/V, (•, ■>*) into (A/ - , (•, ■>*)) such that 

(3.2) [u,v] = (G^u,v)* for all u and v in /C 
(respectively, 

(3.3) [u, v] = (G^lujV)* for all u and v in IC d , 

(3.4) = (G$u,v)* for all u and w in Af). 

Lemma 3.2. The transformations G^\ G^l and G$ are self-adjoint and compact. 
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-<(*) 



Proof. Since [u,v] is symmetric, we immediately get that the transformation G K , G^- 



and G^ are all self-adjoint. For the compactness, we only discuss the case for the trans- 



formation Gfc\ It suffices to show (see, p. 204 of |38|): 
From every sequence {u m } in K(D) which is bounded 

(3.5) ||wm||* < constant, m = 1, 2, 3, • • • , 
we can pick out a subsequence {u m '} such that 

(3.6) (G^\u m > — ui>), (u m i — ui'))* — > when to',/' — > oo. 

Applying Lemma 2.1 with the aid of (13.51) . we find that the sequence {^j 4 -} is precompact 
in L (r e ), so that there is a subsequence {u m >} such that 



d{u m , - Uy 
dv 



2 



ds — > as m , i — >• oo. 



Therefore 

,[{ d(u m , - uy)\ 2 
[u m i - uy , u m ' - uy\ = J Q ^ — j as y as to , / y oo, 

which implies (|3.6p . This proves the compactness of G^ . □ 

Except for the transformations G^\ G^l and Gj$ , we need introduce corresponding 
transformations Gk, G^d and GaA by the inner product (•, •). Since 

(3.7) < [u,u] < (u,u) for all u in K{D) or K d (D) or iV(D), 

there is a bounded linear self-adjoint transformation Gk: of (/C, (•, •)) (respectively, G^d of 
(/C d , (•, •)), Gaa of (A/", (•, •))) such that 

(3.8) [u, v] = (G/cm, v) for all u and u in K 
(respectively, 

(3.9) [it, u] = (G K dU,w) for all u and v in /C d , 

(3.10) [u,v] — (Gj^u,v) for all u and v in A/"). 

Lemma 3.3. The transformations G/c, G^d and G_\f are positive and compact. 

Proof. From [u, u] > for any u S K or /C d or A/", we immediately know that Gic, G^d 
and Gaa are positive. The proof of the compactness is completely similar to that of Lemma 
3.2. □ 

It follows from Lemma 3.3 that Gfc (respectively, G^d, Gjv) has only non-negative eigen- 
values and that the positive eigenvalues form an enumerable sequence {yjc} (respectively, 
{filed}, {/J-jKf}) with as the only limit point. 

Theorem 3.4. The transformations G$ and G/c (respectively, G^l andG^d, G$ and 
Gj\r) have the same eigen] 'unctions. If fi^ and /i/c (respectively, (j,t- d andy^d, yi^j- andyj^f) 
are eigenvalues corresponding to the same eigenfunction we have 



(3.11) 
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( respectively, 

(3.12) M/c d 



1 



H K d 



(3-13) M = j%r)- 

Proof. We only prove the case for the Ga: (a similar argument will work for G^d and 
GV)- Since G^ is positive, we can easily conclude that the inverse (1 + G^)~ 1 exists and 
is a bounded self-adjoint transformation. By virtue of (I3.2|) . (|3.8[) and (|3.ip . we have 



(3.14) = [u > «] = (G JC «,t;) 

= (Gf Jc u > i;)* + (Gx;GWti > «)*, («,«GAC). 

It follows that 

(3.15) Gx: = GW(l + GW)- 1 , 
from which the desired result follows immediately. □ 

Proposition 3.5. Let u and v be two eigen] 'unctions in (K., (•, •)) (respectively, (IC d , (•, ■)), 
(TV, (•, ■))) of the transformation G/c (respectively, G^d, Gj^) at least one of which corre- 
sponds to a non-vanishing eigenvalue. Then u and v are orthogonal if and only if the 4^ I „ 

OV I 1 g 

and |f | r are orthogonal in L 2 e (T g ), that is, 

f du dv , 

3.16 [u,v] = / g— — ds = Q. 

* x e 

Proof. Without loss of generality, we suppose that u is the eigenfunction corresponding 
to the eigenvalue n ^ 0. Then 

[u, v] = (Gku, v) = n{u, v), 

which implies the desired result. □ 

We can now prove 

Theorem 3.6. Let D C (M,g) be a bounded domain with piecewise smooth boundary 
r. Assume that Too is an (n — 1)- dimensional surface inT — T g . If u is an eigenfunction 

of the transformations G^ or Gj$ with eigenvalue fx* ^ 0, then u has derivatives of any 
order in D and is such that 



(3.17) 



( A 2 g u = in D, 
u = on T, 

ga = on T Q0 , A g u = on T - (T e U r 00 ), 
A g u + 7 q |f = on T p , with 7 = j^. 



Proof. Let {uj} be a sequence of functions in K{D) such that \\uj — u\\* — > as j — »• 00. 
We first claim that 

(3.18) u 3 -> u in L 2 (.D). 
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In fact, since Uj £ H$(D) D H 2 (D), it follows from {2/7]) and ([23]) that 



(3.19) 



(3.20) 



/ \uj — ui\ 2 dR < j |V g (uj — ui)\ 2 dR for any j and I, 

Jd Jd 

j \V g (uj - ui)\ 2 dR < ] j \A g (uj - Ui)\ 2 dR for any j and I, 
Jd ^-iW) Jd 



where Xi(D) and A^D) are the first Dirichlet and buckling eigenvalues for D, respectively. 
Since J D \A g (uj - ui)\ 2 dR -> as j, I -> +oo, we find by (|3"H?]) and (pHO]) that |^ - 
ui\ 2 dR — »■ as j,l —> +oo. Therefore the claim is proved. 

For any point p in D, let [/ be a coordinate neighborhood of p, and let E 9 p be a 
bounded domain with smooth boundary such that E C U D -D. Let / be a function in 
Cq{E). Then, by Green's formula (see, for example, p. 6 of [4]), we have 



(u n ,f)*= I (A g u n )(A g f)dR = I u n (A 2 f)dR, 

so that 

(3.21) <«,/>*= / <A 2 J)dR. 



Now by assumption, w = /i*it with /i* ^ and ^ = on T, and hence we have 

(3.22) /)* = /)* = [u, /] = 0. 
Since p is arbitrary in D, it follows from (|3.21[) and (|3.22[) that 

J u{A 2 g f)dR = Q for all / £ Cq(D). 

By applying Green's formula again, we get 

(3.23) / (A g u)(A g f) dR = for all / £ C 4 (L>), 

i.e., u is a weak solution of A 2 u — in D (see [13]). It follows from the interior regularity 
of elliptic equations that u £ C°°(D), and in the classic sense 

(3.24) A 2 g u = in D. 

In exactly the same way, the corresponding result can be proved for G$ . 

Next, suppose that g is continuous. That the boundary conditions of (|3.17[) hold follows 
from Lemma 2.1 and Green's formula. In fact, if 

then uL = and S^L =0, and that 

I r dv I r 00 ' 

/ Q^^-ds = n* ( (A q u)(Av)dR for all v £ K(D). 
J T dv dv J D 



By this and Green's formula (see, p. 114-120 of [23], [2S] and [TU]), we obtain that 

\ j eP^ds= [ (A 2 gu)vdR- [ (A g u)^ds-> r 9{A * U) 
M Jr n ov dv J D J r dv 
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for all v G K{D), where e H' 3 / 2 {T) (see 0). Thus 



(3.25) / (A 2 g u)vdR 



(3.26) 



r-(r e ur 00 ) 



. 1 du\ dv 

(A g u)—ds+ / \ 9 y v ds = 
a;/ J r dv 



for all u G /C(L>). Note that u| = and 



5u I 



= 0, and that 



dv I 



and 



9v I 

dv lr-(r e ur 00 ) 



run throughout space L 2 (T g ) and L 2 (r — (r e UToo)), respectively, when v runs throughout 
space K(D). This implies that 

A g u = on r-(r e ur 00 ), and A g u+— g— = on T e . 

Therefore, (|3.17p holds. In a similar way, we can prove the desired result for G_\f. □ 

Theorem 3.7. Let (M,g) be a real analytic Riemannian manifold, and let D c (M,g) 
be a bounded domain with piecewise smooth boundary T. Assume that T 00 is a real analytic 
(n — 1)- dimensional surface in T — T g . If u is an eigenfunction of the transformations G^l 
with eigenvalue ji* ^ 0, then u has derivatives of any order in D and is such that 



(3.27) 



( A 2 g u = in D, 
u = on T g , 
u = |M = o on r 00 , 

A g u = and = on r-(r e ur 00 ), 

A„u + k p Or- = on r„, with k = -K?. 



Proof. If G Kd u = /i*u, then we have that u = 0on T e and m = ^ = on Too, and that 
(3.28) [ g^^-ds = fi* [ (A q u)(A q v)dR for all v G fc d (-D). 

Applying Green's formula on the right-hand side of (|3.28p . we get that 
,2,,.. JD , f d(A g u) 



(3.29) 



(A z u)vdR- 



di 



v ds 



r-(r e ur 00 ) 



( A 9 u >Q^ds 



/ A„m + — g— — ds = for all u e K (D). 



By taking all u G Cq°(D), we obtain A^u = in D. Note that v L = and u 



Of I 



0, and that v 



and 



dv I 



^lr 00 _ "lr-(r e ur 00 ) — 1 auir-roo ™n throughout the spaces L 2 (r-(r e ur 00 )) 

and L 2 (r — Too), respectively, when i> runs throughout the space K d (D). Thus we have 



A g u = and = on T - (T e U r 00 ), 



Am 



on I\ 



□ 
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Theorem 3.8. Let (M,g), D and Too be as in Theorem 3.7. Assume that and Kk are 
the k-th Steklov eigenvalues of the following problems: 



(3.30) 



in D, 



A 2 g u = 
u = on T g , 
u=%± = on r 00 , 

|H = o and =0 onT-(T e D T 00 ), 

?f>|£=0 on r„ 



AgU 



and 



(3.31) 



u = on T 



m D, 



on r 00 , 



A g u = and 
Agu + Kg^ = on T g , 



o on r - (r e u r 



respectively. Then ^ < K k for all k > 1. 

Proof. For < a < 1, let Uk = Uk(a,x) be the normalized eigenfunction corresponding 
to the fc-th Steklov eigenvalue A& for the following problem: 



(3.32) 



A^gUk 



in D, 



Uk 

aA g Uk 



on T c 







dui. 



(1 - a 



on r 00 

dv 



= and 



= on r - (T g u r 



(II) ) 



Ag Uk + A q ^ = on r 



It is easy to verify (cf, p. 410 or Theorem 9 of p. 419 in [6]) that the fc-th Steklov eigenvalue 
A/c = A/ £ (a) is continuous on the closed interval [0, 1] and differentiable in the open interval 
(0, 1), and that u k (a, x) is also differentiable with respect to a in (0, 1) (see, [H]). We will 
denote by ' the derivative with respect to a. Then 



(3.33) < 







u' k = on Tg, 

— = on 



in D, 



u k = 
A„u k 



n,. 



00, 



a 



A 



9 U k 
du k 



^ + (l-a)^ =0and 



dui 



d(A g 



I AX + A'^+A?|=0 on r B . 



= on r - (r e u r 00 ) 
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Multiplying p.33|) by Uk, integrating the product over D, and then applying Green's for- 
mula, we get 







(A 2 g u' k )u k dR= f (A 2 g u k )u' k dR- f (A g u k )^ds 

ds+f (A g u' k )^±ds 

JdD av 



, d(A g u k ) 
u k ^ ds - 

8D dv 



, a % du', , 

(AgU k )— ±dS 

av 



„ 9(A g u' k ) 
Uk ^ — 

dD OV 



r-(r 8 ur 00 ) 



du k \ du 



dv 







' k ds 



r-(r e ur 00 ) 



r-(r e ur 00 ) 



{A g u k )—±ds 
ov 

\A g u k )—ds 



-Xq—^--\q—^- ^-ds 
1 ' dv 



1 — a. du k \ du' k 
a dv J dv 



K i n du k x n du' k \ du k 



-A' 
-A' 



r-(r e ur 00 ) 
du k 
dv 



AgU k , 

a a dv 



1 du k 1 — a du' k \ du k 



o 



ds 



( duk_ 
\ dv 



r-(r e ur 00 ) 



r-(r e ur 00 ) 



a dv J dv 
1 — a \ du k 1 du k 



ar 

1^ du k 

a dv 



d 



v a dv 



dm 



ds 



ds, 



i.e. 



A' fc («) = 



Jr-(r e ur 00 ) (I dv) ds 
Ir^m'ds 



> for all < a < 1. 



This implies that X k is increasing with respect to a in (0, 1). Note that if we change the a 
from to 1, each individual Steklov eigenvalue X k increase monotonically form the value 
Ob which is the fe-th Steklov eigenvalue of p.30[) to the value n k which is the k-th Steklov 
eigenvalue (|3 .3 1[) . Thus, we have that g. < n k for all k. □ 



Conversely, the following proposition shows that a sufficiently smooth function satisfying 
(|3.17[) (respectively, (|3.27[1 ) is an eigenfunction of G$ or G$ (respectively, G^l). 

Proposition 3.9. Let D be bounded domain with piecewise smooth boundary. Assume 
that u belongs to C A {D). 

a) IfTg^T and u satisfies \3. 17\ ), then u G /C and u is an eigenfunction of G^ with 
the eigenvalue fi* = 7 _1 , 



(3.34) 



G ( * ] u = 7 ~V 



b) IfTg^T and u satisfies \3.27\ ), then u G K d and u is an eigenfunction of Gj^l with 
the eigenvalue /i* = k^ 1 , 



(3.35) 
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c) IfV e = T and u satisfies then u € M and u is an eigenfunction of G$ with 

the eigenvalue /Lt* = 7 , 

(3.36) G% ] u = 7~V 

Proof, i) T e ^ r. We claim that there is no eigenvalue 7 = 0. Suppose by contradiction 
that there is a function u in C (D) satisfying 

( A 2 g u = in D, «=0on T, 
[ ' \ % = on r 00 , and A g u = on T - r 00 . 

By multiplying the above equation by u, integrating the result over D, and using Green's 
formula, we derive 

0= f u{A 2 a u)dR= [ \A„u\ 2 dR- [ u d ^ u ) ds 
Jd Jd Jt ov 

+ [ (A g u)pds = [ \A g u\ 2 dR. 

JT V v JD 

This implies that A g u = in D. Since u = on T, by the maximum principle we get that 
u — in D. The claim is proved. 

In view of assumptions, we see that u € K.. By (|3.17p and Green's formula, it follows 
that for an arbitrary v G K(D) 



{G y K 'u,v)* = [u,v]= Qjr-j-ds 



-7- 1 / {A g u)^ds=- 1 - 1 [ (A g u)^ds 



1 ■ dv J r dv 



-7- 1 



vds ~ I (& g u)(&gv)dR+ I v{A 2 g u)dR 



IT vi/ Jd JD 

7" 1 / (A g u)(A g v)dR = 1 - 1 {u,v)\ 



D 



Therefore, 



(G% ] u - j^u, v)* = for all v e K(D), 



which implies (|3.34l) . By a similar way, we can prove b). 

ii) T e = T. We claim that there is no eigenvalue 7 = 0. If it is not this case, then there 
is a function u in C 4 (D) satisfying 

A 2 g u = in D, 



Setting v := A g u in D, we get 



u = on r, 
A g u = on T. 



A g v = in D, 



v = on T. 

By the maximum principle it follows that v = in D. Thus, we have 

A g u — in D, 



u = on r, 



so that u — in D. 
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Now, if u is a solution of p,17[) with eigenvalue 7 > 0, proceeding as in a), we can prove 
that u G TV and ([336]) holds. □ 

Remark 3.10. Each of transformations G^, G* Kd and Gjy corresponds to a biharmonic 
Steklov problem given by the quadratic forms 



(u,u)* = / \A g u\^dR 
Jd 

and 

9u x 2 



M] = /r. e W ^ 

and the function classes of /C*, /C rf * and A/"*, respectively. The eigenvalues 7^ and Kfc of 
these biharmonic Steklov problems are given by 

(3.38) 7fe and Kk — 1//4 k = 1,2,3, 

Since is the only limit point of /j,t, the only possible limit points of jk and Kk are +00. 



4. Biharmonic Steklov eigenvalues on an ?i-dimensional rectangular 

parallelepiped 

Let D = {x e K n |0 < Xi < U, i = 1, • • • , n} with boundary T, and let L+ = {x G 
R n |0 < < ij when i < n, x n = 0}. Let T ln = {x G K"|0 < Xi < h when z < n, x„ = /„}. 
Our first purpose, in this section, is to discuss the biharmonic Steklov eigenvalue problem 
on n-dimensional rectangular parallelepiped D: 

A 2 u = in D, 
U = on r, 



(4.1) 



|H = on I 11 ", Au = on V - (L+ U L z "), 



Am + jg^j = on T+ g = constant > on L+. 
We first consider the special solution of (|4.1|) which has the following form: 

u = X(xi, ■ ■ ■ ,x n -i) Y(x n ). 

Since 

Au = (A„_iX(a;i,--- ,x n -i))Y(x n ) + 2VX(xi,- ■ ■ ,s„_i)- VF(a; n ) 
+ ,ss n _i))y"(aj n ) = (A„_ 1 X(a; 1 , • • • ,ar n _ 1 ))y(x„) 

+ (A>i,-- - .Sn-ijjy"^) 

and 

A 2 u = (A^ l _ 1 X(xi, ■■■ , x n -i))Y(x n ) + 2(A„_i • • • , x n -ij)Y" {x n ) 
+ (X(x u --- ,x n ^))Y""{x n ), 



where 



A»_iX(a;i, • • ■ ,x n -i) = 



1 
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we find by A 2 u — that 

X\ , * • • , Xn— 

i))Y(x n ) + 2(A„_iX(xi, • • • , x n -i))Y"(x n ) 
+ (X(x u --- ,x n - 1 ))Y""{x n )=0, 

so that 



(4.2) 



A^Xixi,--- ,x n -i) A n ^ 1 X(x 1 , • • • ,x„_i) Y"(x n ) , F""(xn) 



X(xi, • • • ,x„_i) " X(xi, ■ ■ ■ ,x„_i) y(x„) Y(x„) 
Differentiating (|4.2I) with respect to x„, we obtain that 



0. 



2 A n _iX(xi, ■ ■ ■ ,x n _i) 
X(xi,- • • ,x„_i) 
The above equation holds if and only if 



Y"(x n ) 



Y(x n ) 



Y""{x n ) 



Y(x n ) 



= 0. 



(4.3) 



A n _ 1 X(x 1 , ■ ■ ■ ,x n _i) 
X(xi, ■ ■ ■ ,x„_i) 



Y""(x n ) 
Y(x n ) 



Y"(x n ) 



-7] 



where r] 2 is a constant. Therefore, we have that 

(4.4) A n _iX(x u --- ,x n -i) +r 1 2 X{x 1 ,--- ,x n _i) = 

and 



Y""(x n ) 



Y(x n ) 



2r? 



Y"{x n ) 



Y(x n ) 



= 0. 



From (|4.4p . we get 

(4.5) A^_ 1 X(xi,-- - ,x„_i) = -i 1 2 A n - 1 X{xi,--- ,x n -i) = ?? 4 X(xi,- 
Substituting this in (|4.2j) . we obtain the following equation 

(4.6) y""(x n ) - 2^ 2 r"(x„) + rfY{x n ) = 0. 

It is easy to verify that the general solutions of (|4.6I) have the form: 

(4.7) Y(xn) = A cosh rjx n + B sinh rjx n + Cx n cosh r\x n + Dx n sinh rjx 
By setting Y(0) = y(Z n ) = 0, F'(0) = 1, Y'(/ n ) = 0, we get 

— 77Z 2 ^ ^ sinh 2 77? 



, x„_i). 



(4.8) F(x n ) 



sinh z 7yZ„ — rfl 2 



sinh ?7X r , 



?7^n — (sinh r]l n ) cosh r?i T . 
sinh 2 ?7i„ — ?7 2 Z^ 



sinh r\l n — ?/ 2 Z 2 
x n sinh?7a; n . 



coshryx,, 



It is well-known that for the Dirichlet eigenvalue problem 

, . f A„_iX(xi, • • • ,x„_i) + ?y 2 X(xi, ■ • • ,x„_i) = in f2, 

\ u = on <9{(xi, • • • ,x n _i)|0 <Xi <k, i = 1, • • • ,n - 1}, 

there exist the eigcnfunctions 



(4.10) A(xi, • • • ,x n _ij = c I sm — — xi 

which correspond to the eigenvalues 



. m„_i7r 
sm — x„_i 

In— 1 



2 = 1 



1 ~T~ 



where rrii = 1, 2, 3, • • • 
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Therefore, 
(4.11) u 



(x(x 1 , ■ ■ ■ ,x n _i))y(x„) 



rriiTT 
c | sin — ; — x\ 



h 

sinh 2 nl n 



sinh rjl n — i] 2 l 2 
ilUi — (sinh rjl n ) cosh rjl n 



m n _i7r 
sin — x n ~i 

COSh rjX n 



-nil 



sinh nl n — rfl n 



sinh nx ri 



x n sinh nx n 



Since 



we obtain 



F"(0) = 2rj 



sinh r^„ — n 2 l 2 

V^-n ~ (sinh ryZ„) cosh i]l n 



sinh 2 r^„ — rj 2 l n 



and F'(0) = 1, 



i„=0 



(A n _iX(a;i, ■ • • ,x n - 1 ))Y{0) + (X(x l ,'-- ,a;„_ 1 ))Y"(Q) 
?7^i — (sinh cosh 



2// 



and 
so that 

with 



du 
dv 



7 



sinh 2 7]l n — rj 2 l n 
= X(xi, ■ ■ ■ ,x n -t), 

+ 
p 

du 

Ati + = on T+ 

ov e 

2r)l n ( (sinh?7/ n ) coshr,^ - rjl n 
gin \ sinh 2 r)l n - ifll 



X(x\, • • • , x n -i), 



Our second purpose is to discuss the biharmonic Steklov eigenvalue problem on the 
n-dimensional rectangular parallelepiped D: 

A 2 u = in D, 



(4.12) 



u = on r+, u = |f = on r ; » 

|2± = ^2 = on r - (r+ u rM, 



Am + S£?|f = on T+, g = constant > on T+. 



Similarly, (14.121) has the special solution u — (X(xi, • ■ • , x n -i)) Z(x n ) with Z(x n ) having 
form (|4.7p . According to the boundary conditions ol (I4.12j) . we get that the problem (|4.12l) 
has the solutions 



u(x) = u{x u ■■■ ,x n ) 



TO17T 
C I COS — X\ 



TUn-lTT , , . 

COS — X n -i Z[x n ), 



h J \ ln-1 

where mi, • ■ ■ , m n -i are whole numbers, and Z(x n ) is given by 



(4.13) Z(x n ) = 



sinh^ j3l n - /3 2 l n 



sinh f3x n 



sinh pi n 



sinh^ pi n -P 2 11 



X n COSh Px n 



pi n — (sinh pi n ) cosh /3Z„ 
sinh 2 /3Z„ - p 2 l 2 n 



x n sinh 
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P= [E^CnWW 8 ] V2 with ET^mi ^0. Since ^| r + = ••■ , x^), (Au))^ = 

(*(*!, • • • , x n -i))Z"(0) and Z"(0) = 2ry ( ^"^K^^'" ) , we get A U + = on 
r+, where 

_ 20ln / (sink ffl n ) cosh ffl n 
? ~ V sinh 2 /?/„ - /3 2 Z2 J" 



5. Asymptotic distribution of eigenvalues on special domains 

5.1. Counting function A(r). 

In order to obtain our asymptotic formula, it is an effective way to investigate the 
distribution of the eigenvalues of the transformation Gx (respectively, G^d, Gj^f) instead 
of the transformations G^ (respectively, Gj^l, G$). It follows from p. Ill) — p,13[) and 
(|3.38|) we obtain 

(5.1) fJL k = (l + X k )-\ k= 1,2,3,- •• , 

where /x& denote the fc-th eigenvalue of Gk, or G^d or Gaa, and j- is the fc-th eigenvalue of 
Gfc or G^l or G$) (More precisely, = 7^ for Gj£ and G$, and Afc = Kk for G^]).) 
Since -A(r) = J2\ k < T we nave 

(5.2) A(r)= J] L 

5.2. D is an n-dimensional rectangular parallelepiped and gfy = 5^. 

Let -D be an n-dimensional rectangular parallelepiped, gij — &y in the whole of D, 
g = constant > on one face r+ of the rectangular parallelepiped, i.e., D = {x <E 
R"|0 < Xi < h, i = 1, ••• ,n}, T+ = {.t e M"|0 < x t < h when z < n, x„ = 0}) and 
Too = r'" = {x 6 M"|0 < < Z,; when i < n, x n = l n }. Without loss of generality, we 
assume l{ < l n for all i < n. 

For the above domain D, except for the K(D) and K d (D) in Section 3, we introduce 
the linear space of functions 

— (JU 

K°(D) = {u\u e Hq(D) n H 2 (D) n C°°(D), — = on r 00 , Au = on T - (T+ U r 00 )}. 
Clearly, 

(5.3) K°(D) C K(D) c K d (D). 

Closing K°, K and respect to the norm ||tt|| = y/ (u, u), we obtain the Hilbert spaces 
JC°, K and lC d , and 

(5.4) K° C K C /C d . 
According to Theorem 3.3, we see that the bilinear functional 

, f du dv . 
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defines self- adjoint, completely continuous transformations G°, G and G d on KP, K, and 
K d , respectively (cf. Section 3). Obviously, 

(G°u,v) = (Gu,d) for all u,v in K°, 

{Gu, v) = {G d u, v) for all u, v in K, 
from which we deduce immediately by Proposition 2.3 that 

(5.6) M°<Mfc</4, k = 1,2,3,- ■■ , 

where and {^f} are the eigenvalues of G° and G d , respectively. Hence 

(5.7) A°{t) < A(t) < A d {r) for all r, 
where 

(5.8) A°(r) = J2 1 



and 
(5.9) 



A d (r) 



M^>(1+t) 



We shall estimate the asymptotic behavior of A°(t) and A d {r). It is easy to verify (cf. 
Theorems 3.6, 3.7) that the eigenfunctions of the transformations G° and G d , respectively, 
satisfy 

A 2 w = in D, 
u = on r, 

and Au = on r - (T g U r ; »), 
g = constant > on ^ 



(5.10) 



|a = on r z - 

Am + jg?p = on T+, g = constant > on Tj. 



and 



(5.11) 



A 2 u = in D, 



on 



Au + K g^ 



<)u 
dv 



on P™, 



and Au = on T-(r+ur'"), 
Q 

As being verified in Section 4, the functions of form 



constant > on T+. 



(5.12) 



u(x) 



miir 
sin — — xi 
h 



sin ■ 



-17T 



r(Zn) 



are the solutions of the problem (|5 . 10[) . where mi,-- - ,m„_i are positive integers, and 
y(x„) is given by (|4.8[) . Since the functions in (|5.12[) have derivatives of any order in 
D, it follows from Proposition 3.9 and Theorem 3.4 that they are eigenfunctions of the 
transformation G° with eigenvalues (1 + 7) , where 

rn-i 1 1 / 2 

2nl n ( (sinh?7/ n ) cosh?7Z„ — nl n 



(5.13) 



7 



Note that the normal derivatives 
du 
dv 



sinh 2 nl n — n 2 l 2 



V 



n-l 
,i=l 



(5.14) 



mi7r 
sin — — X\ 



h 



_17T 



- X n —\ 



when mi, • • • , m n _i run through all positive integers (see, Section 4), form a complete sys- 
tem of orthogonal functions in L 2 (T+). It follows from Proposition 3.5 that if mi, ■ ■ • , m n _i 
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run through all positive integers, then the functions (|5.12|) form an orthogonal basis of the 
subspace of KP, spanned by the eigenfunctions of G°, corresponding to positive eigenvalues. 
That is, when mi, • ■ ■ , m n _i run through all positive integers, then (1 + 7) _1 , where 7 is 
given by (|5.13[) . runs through all positive eigenvalues of G°. 

Similarly, for the problem (|5.11l) . the eigenfunctions {«/,} of the operator G d on /C d , 
corresponding to non-zero eigenvalues, form an orthogonal basis of the subspace of K d . The 
non-zero eigenvalues of G d are u d = (1 + Kk)^ 1 , where Kk is the fe-th Steklov eigenvalue of 
(|CT|l . 

In order to give the upper bound estimate of A d (r), we further introduce the following 
Steklov eigenvalue problem 

A 2 u = in D, 



(5.15) 



u = on r+, u = = on T l », 
= and = on T - (T+ U r'« 



at/ v e 

Ait + = on T+, g = constant > on T+. 

Let Cfc be the fc-th eigenvalue of (|5.15|) . By Theorem 3.8, we have 

(5.16) < Kfc for all k > 1. 
We define 

(5.17) ^ = TT^' a/(t)= S l 

It follows from ([5~16|) and (j5T7)) that 

(5.18) A d {r) < A f (r) for all r. 

We know (cf. Section 4) that the problem (15.151) has the solutions of form 

(5.19) u{x) = c ^cos -j— x ^j ' • ' ^ cos ~y ~ x n-ij Z{x n ), 

where mi,-- - ,m„_i are non-negative integers with J^ILi m * 7^ 0' ari< ^ Z{x n ) is given 
by (14. 13[) . This implies that if mi, - - • , m n _i run through all non- negative integers with 

HlZi m i 0) tncn 

i 1/2 



_ 2ffl» / (sinh ffi») cosh - /3l n . 



sinh z 0i n - pHl 



n-1 

.2 — 1 



runs throughout all eigenvalues of problem (|5 . 1 5[) . 

We first compute the asymptotic behavior of A^(j). By (I5.17p . (|5.20[) and the argument 
as in p. 44 of [50] or p. 373 of [6] or p. 51-53 of [39], -<4^(t) =the number of (n — l)-tuples 
(mi, • • • ,m„_i) satisfying the inequality 



2ffl n / (sinh /3l n ) cosh ffl n - f3l n \ 
gin V sinh 2 pi n - pHl ) 



(5.21) v , 2:7 j <^ 

where mi, ■ ■ • , m n _i are non-negative integers with ^"Ji 1 m * 7^ 0- By setting 

„ , / (sinh s) cosh s 

(5.22) t(«) = 2s y —^- 2 

\ smn s — s z 

we see that 

lim t(s)/s = i 

s— >-j-oo 
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We claim that for all s > 1, 

, [— 3s(sinh 2 s) + 3s 2 (sinh s) coshs + (sinh 3 s) coshs — s 3 (sinh 2 s + cosh 2 s) 

t(s) = 2 



In fact, let 



Then 



(sinh 2 s — s 2 ) 2 

-3s(sinh 2 s) + 3s 3 (sinh s) cosh s + (sinh 3 s) coshs — s 3 (sinh 2 s + cosh 2 s) 
9(1) > 0, and 



> 0. 



9'(s) = — 3(sinh 2 s) — 4s 3 (sinh s) coshs + 3(sinh 2 s) cosh 2 s + sinh 4 s 
= 4(sinh s) [sinh 3 s — s 3 cosh s] > for s > 1 , 

This implies that 9(s) > for s > 1. Thus, the function t(s) is increasing in [l,+oo). 
Denote by s = h(t) the inverse of function t(s) for s > 1. Then 

h(t) _ 1 



lim 

t->+oo t 2 



Furthermore, we can easily check that 



(5.23) 



Note that, for s > 1, the inequalities t(s) < t is equivalent to s < ft,(t). Hence (|5.2ip is 
equivalent to 



which can be written as 

n-l 

(5.24) X)KA) 2 < 



7rZ„ 



< h(gl n r), 



h(gl n r) 



rrn = 0,1,2,- 



We consider the (n — l)-dimensional ellipsoid 

n-l 



7T<„ 

Since A* (r) + 1 just is the number of those (n — l)-dimensional unit cubes of the z-space 
that have corners whose coordinates are non-negative integers in the ellipsoid (see, VI. 
§4 of [6]). Hence A'(t) + 1 is the sum of the volumes of these cubes. Let V(t) denote 
the volume and T(r) the area of the part of the ellipsoid situated in the positive octant 
Zi > 0, i = 1, • ■ • , n — 1. Then 

(5.25) V(r) < A f {r) + 1 < V{t) + (n - l)*T(r), 

where (n — 1)^ is the diagonal length of the unit cube (see, [6] or [39]). Since 



{n-l) 



by f[5723]) . we get that 

(5.26) V(t) - w n _i(47r)-("- 1 )/ 1 ---/„_ l£ .™- 1 r"- 1 +0(r™- 2 ), as r -> +oo. 
Note that 
(5.27) 



T(t) ~ constant • r 



n-2 
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It follows that 

lim =uj n ^ 1 {An)-^- i n i --4 n - 1 Q n -\ 

i.e., 

(5.28) Af{r) ~ ir+l^-V"- 1 , as r +00, 

where |T+| denotes the area of the face T+. 
Next, we consider A°(t). Similarly, 

2^n ( (sinh?7^„) cosh?7Z„ — r]l T: 
gl n \ sinh 2 r\l n - jflj. 



(5.29) ^ v ' 7 ' " < r, 



is equivalent to 

Tjl n < h(gl n r), 

i.e., 



£[^A] 2 <(^V mi = i, 2)3 

i— 1 ^ n / 



Similar to the argument for A-^(r), we find (see also, [25] or §4 of [6]) that 



n-l 



, \ ^ ,mj.2 ( h(gl n T) ^2 { 
#{(%."• ,m n -i)\ 2^ H-) <( ; ) , m, = 1,2,3, ••• ) 

2=1 

k-'n — 1 i-pi. 1 n— 1 t),— 1 



(5.30) lim = -^fi-lr+le"" 1 . 

v ; t"" 1 (47r)("- 1 ) e 

Noting that q = on T e - T+, by (^7f| . ([5TTB]) . (f^28|) and ((ODjl . we have 

(5.31) A(t) ~ / e n-1 da asT4 +00. 



5.3. A cylinder I? whose base is an n-polyhedron of M™ 1 having n-l orthogonal 
plane surfaces and g. L j = o~,-j. 

Lemma 5.1. Let = x [0,Z„], r = 1,2, w/iere = {(a?!, • • • ,x n -i) G R™ _1 | a?i > 
/or 1 < i < 71 - 1, and X)"=i ff - 1 }' arid r e 

is an (tj — 1) -dimensional cube with side 
length I — maxi<i<„_i l{. Assume that Iqq = Tg x {l n }, r = 1, 2. Assume also that g is 

(r) 

a positive constant on Ty, r = 1,2. I/Z < Z n , iften 

(5.32) «£p (1) )>«£(£> (2) ) /or A; = 1,2,3, • • • , 

where <;l(D^) (similar to c 0/ 15*. m Theorem 3.8) is the k-th Steklov eigenvalue for 
the domain . 
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Proof. Let be the &;-th Neumann eigenfunction corresponding to oi£ ' for the (n — 1) 



(r) 



(r) 

dimensional domain T g , (r = 1,2 



(5.33) 
Put 



(r) 



i.e., 

M.,(r) _ 



A' 



(r) 



in T 
on dr 



(r) 

e ■ 



,(<•) 



(ar) = («W(si,.-. ,a: n _i))(zW(a; n )) in £>« 



where Z^ r '(x n ) is as in (|4. 13|) with /3 being replaced by 
r (x) satisfies 

in D {r \ 



It is easy to verify that 



(5.34) 



with 



A' 2 



= 



4 r) = o 



on r 
on r 



(r) 

e i 

(r) 



A^+^( J DW)^=0 



00 ' 

on at>M 
on r". 



(5.35) 



(sinh 



'a^ 2 n )cOsh \/Q!fe i« 



sinh 



/ 2 



(1) 



It follows from p. 437-438 of [6 that the fc-th Neumann eigenvalue a 
is at least as large as the fc-th Neumann eigenvalue a k for the domain T e 



for the domain L^ 1 "* 
Recalling that 



2s 



(sinh s) cosh s — s 
sinh 2 s — s 2 



is increasing when s > 1, we get 

&DW)>4(DW), fc = 1,2,3,.-. 



if 2 < 2 n . Here we have used the fact that 



> 1 since any Neumann eigenvalue 

for has the form Yn=i (^^p")^' ^ n °ther words, if I < l n , then the number A? (t) of 
eigenvalues less than or equal to a given bound r for the domain is at most equal to 
the corresponding number of eigenvalues for the domain D^ 2 \ □ 

Similarly, we can easily verify that the number A ' (t) of eigenvalues less than or equal 
to a given bound r for an arbitrary n-dimensional rectangular parallelepiped D is never 
larger than the corresponding number for an n-dimensional rectangular parallelepiped of 
the same height with its base an (n — l)-dimensional cube whose side length is at least 
equal to the largest side length of the base of D. 



5.4. D is a cylinder and = dy. 

Let D be an open n-dimensional cylinder in R™, whose boundary consists of an (n — 
l)-dimensional cylindrical surface and two parallel plane surfaces perpendicular to the 
cylindrical surface. Assume that gij = 5ij in the whole of D, that T g includes at least 
one of the plane surfaces, which we call Tj, and that g is positive constant on and 
vanishes on T e — T+. We let the plane surface be situated in the plane x n = and let 
another parallel surface T ln be situated in the plane {x € M" \x n = l n }. We now divide 
the plane x n = into a net of (n — l)-dimensional cubes, whose faces are parallel to the 
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coordinate-planes in i„ = 0. Let Ti, • • • ,T P be those open cubes in the net, closure of 
which are entirely contained in Tj, and let Q p +i, • • • , Q q be the remaining open cubes, 
whose closure intersect Tj. We may let the subdivision into cubes be so fine that, for every 
piece of the boundary of T~t which is contained in one of the closure cubes, the direction 
of the normal varies by less than a given angle z9, whose size will be determined later. 
(This can be accomplished by repeated halving of the side of cube.) We can make the 
side length I of each cube be less than l n . Furthermore, let Dj, (j = 1, • • • , p), be the 
open n-dimensional rectangular parallelepiped with the cube Tj as a base and otherwise 
bounded by the "upper" plane surface r ; ™ of the cylinder D and planes parallel to the 
coordinate-planes x\ = 0, • • • , £„-i = (cf. [39]). 
We define the linear spaces of functions 

K = {u\u G Hl{D) n H 2 {D), ^ = on r 00 }, 

= { Uj \ Uj G H^Dj) n H 2 (D 3 ) nr^),!^ on rj», and 

Au 3 = on dD j - (T, U ty)}, (j = 1, • • • ,p) 

with the inner products 

(u, v) = J D (Au)(Av)dR + f Tg q% f ds for u, v G K, 

{uj,Vj)j = (uj^j)* + [u^v^j = J D . (Auj)(Avj)dR + J r , ^ds for u 3 ,v 3 e K$, 

respectively. Closing K and K® with respect to the norms \\u\\ = y/ (u, u) and = 
yj (uj,Uj)j, we obtain the Hilbert spaces /C and K.® (j = 1, ■ • • ,p), respectively. Clearly, 
the bilinear functional 

f du dv 

[u,v\= Qjr~A~ ds 

Jr e dv dv 

= J r Q ds > (i = 1» • • • ,P), 

define self-adjoint, completely continuous transformations G and on /C and K,® by 

(5.36) (Gu,v) — [u,v] for u,v in /C, 

(5.37) (G < jiij,Vj) j — [uj,Vj]j for Uj and Vj in /C°, 

respectively. By defining a space 
p 

K? = J2 ® JC °3 ={u°\u° =m + --- + u p , Uj G K*} 

3=1 

with its inner product 

p 

(5.38) (u°,v°)=J2(u j ,v j ) j , 

3=1 

we find that the space /C° becomes a Hilbert space. If we define the transformation G° on 
/C°by 

(5.39) G°u° = G° 1 u 1 + --- + G° p u p for u° = u x + ■ ■ ■ + u p in K°, 

we see that G° is a self-adjoint, completely continuous transformation on /C°. If we put 

(5.40) [«°,« ]=f>i,«i]i> 

3 = 1 
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wc find by (IQ71) — (jOO]) that 

(5.41) {G°u°,v ) = [u°,v°] for all u° and v° in K° . 

Let us define a mapping of /C° into /C. Let u = u\ + ■ • • + u p , uj S -ffj*, be an element of 
/C° and define 

(5.42) u = n°u°, 

where u{x) — Uj(x), when x G Dj, and u(x) — 0, when x E D — L)^ =1 Dj. Clearly u G /C 
and thus (|5.42|) defines a transformation II of Kf{ © • • • © /C° into /C. It is readily seen that 

(5.43) [IlV,nV] = [u°,v°] for all u° and w° in K. . 
and 

(5.44) (G%°,v ) = (Gn% ,n% ) for all u° and v° in /C°. 
By (I5.43|) and (|5.44l) . we find by applying Proposition 2.3 that 

Mfe < for fc = 1,2,3,- •• . 

Therefore 

(5.45) A°(r) < A(t). 
The definition of G° implies that 

(5.46) G°/C°c/C°, (7 = l,-.-,p), 
and 

(5.47) G°m° = G°u°, when w° g /C?. 
From (t5T^D|) . (jiplTI) . (flOS)) . (1CT7)) and Proposition 2.4, we obtain 

(5.48) A°(t)=£^(t), 

where A° (r) is the number of eigenvalues of the transformation G° on JC® which are greater 
or equal to (1 + t) . Because Dj, (J = 1, • • • ,p), is an n-dimensional rectangular paral- 
lelepiped we find by (|5.30l) that 

(5.49) A°(t) ~ w n _ 1 (47r)-("- 1 )|r j | £ -"- 1 r™- 1 as r ^ +oo, 

where \Tj\ denotes the area of the face Tj of Dj. By (|5.48[) and (|5.49|) we infer that 

p 

(5.50) A°(t) ~ w„_ 1 (4^)- ( "- 1) ^ ir.l^-V™- 1 as r ^ +oo. 

i=i 

Next, we shall calculate the upper estimate of A{t). Let Pj, (j = p + 1, ■ ■ ■ , q), be the n- 
dimcnsional rectangular parallelepiped with the cube Qj as a base and otherwise bounded 
by the "upper" plane surface T ln of the cylinder D and planes parallel to the coordinate- 
planes x\ = 0, • • • , x n ^i — 0. The intersection Pj O D is a cylinder Dj,(j = p + 1, • • • , q), 
with f j := Qj (~l f + as a base. Obviously 

(5.51) D = J2Dj. 
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We first define the linear spaces of functions 

K d = {u\u G H 2 (D), u = on T e , u = = on r'™)}, 
Kf = { Uj \ Uj G H*(Dj), Uj = on r„ uj =^=0on r<"}, (j = 1, • • • ,q) 
with the inner products 

(5.52) («, ») = ^(AtiJCA^dR + / re e^^ <fa, 
and 

(5.53) (u j ,v j ) j = [ (Au 3 )(Av 3 )dR+ [ g^^ds, 

JDj JTj OV OV 

respectively. Closing K d and K d with respect to the norms ||u|| = \J (u, u) and = 

\J (uj,Uj)j, we obtain Hilbert spaces JC d and 1Cj,(j = 1, • • • , <?), and then we define the 
Hilbert space 

Q 

(5.54) /C d = ^©/C J rf = {u d |u d = ui + --- + M g , Uj6X^} 

j'=i 

with its inner product 

(5.55) (u d ,v d )=J2(uj,v j ) j . 
The bilinear functional 

(5.56) [uj,Vj]j = [ g^-^-ds, (j = 1, • • • , <?), 

define a self-adjoint, completely continuous transformation G d on JCj given by 

(5.57) (G d Uj,Vj)j — [uj,Vj]j for all Uj and vj in K d . 

The self-adjoint, completely continuous transformation G d on K. d is defined by 

(5.58) G d u d = Gf Ul + ■■■ + G d u q for u d = Ul + • • • + u„ in /C 1 . 
With 

(5.59) [tt d ,w d ] = Y^ u i-> v i\h 

3 = 1 

we find by (|535j) . (I537D — (153^1 ) that 

(5.60) (G*V,v d ) = for all u d and w d in /C d . 
Let us define a mapping II of JC into AC d . Let u G K(D) 1 and put 

M d = LTu = Ul + • ■ • + u q , 
where Uj[x) — u(x), when x G Dj. It can be easily verified that 

(5.61) {Ilu,ilw} = (u, v) for all u and v in /C. 
and 

(5.62) (Gu,v) = (G d Uu,Uv) for all u and v in /C. 
From (|5.61|) — (|5.62|) . with the aid of Proposition 2.3, we obtain 

Mfe <4 for A = 1,2,3,--- , 
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and hence 

(5.63) A(t) < A d {r). 

By G d K. d clC d , (J = 1, •••,?), and G d u d = G d u d when u d € /C^, we get 

(5.64) A d (T)^4(r), 

3=1 

where A d {r) is the number of eigenvalues of the transformation G d on K, d which are greater 
than or equal to (1 + r) _1 . Further, we define Aj(t) similar to ()5.15j) and (|5.17|) . i.e., 

4(r) = J2 1 with 4 = TTT' 

* 1 T" Sfe 

where Cfc is the fc-th Steklov eigenvalue of the following problem 
A 2 Uj = in Dj, 
u-j = on Tj, ■ = ^ = on 

f^ = T = ° on ^--(r.-ur^), 

Auj + SQ-^j- — on Tj, g — constant > on T+. 
From Theorem 3.8, it follows that 

Qc < for all fc > 1, 

and hence 

(5.65) A d {r) < Aj(r) for all r and j = !,■■■ ,q, 

where ^ is the fc-th eigenvalue of the transformation . Since Dj,(j — 1 , • • • , p) , is an 
n-dimensional rectangular parallelepiped, we find from (|5.28l) that 

(5.66) 4(r)^ Wn _ 1 (4 7 r)-(™- 1 )|r j | e "- 1 r"- 1 , (j = l,---,p). 

It remains to estimate Aj(r),(j > p + 1). According to the argument in p. 438-440 
of IB], each of the (n — l)-dimcnsional domains Tj is bounded either by n — 1 orthogonal 
plane surfaces of the partition (the diameter of the intersection of any two plane surfaces 
lies between I and 3Z), and an (n — 2)-dimensional surface of the boundary (see, in two 
dimensional case, Figure 5 of p. 439 of |BJ), or by 2n — 3 orthogonal plane surfaces of the 
partition (the diameter of the intersection of any two plane surfaces lies between I and 31), 
and a surface of the boundary dT g (see, in two dimensional case, Figure 6 of p. 439 of [B]). 
The number q — p is evidently smaller than a constant C/l n ~ 2 , where C is independent of 
I and depends essentially on the area of the boundary dT g . Now, we take any point on the 
boundary surface of Tj and take the tangent plane through it. This tangent plane together 
with the plane parts of dTj bounds an n-polyhedron of R™ -1 with a vertex at which n — 1 
orthogonal plane surfaces meet (see, Figure 5 of p. 439 of [BJ in two dimensions), e.g., if d is 
sufficiently small it forms an (n — l)-dimensional rt-polyhedron of R™ with a vertex having 
n — 1 orthogonal plane surfaces (the diameter of the intersection of any two plane surfaces 
is also smaller than 4/), or else an (n — l)-dimensional 2(n — l)-polyhedron of R™ -1 (see, 
Figure 6 of p. 439 of [BJ in two dimensional case), the diameter of the intersection of any two 
plane surfaces (except for the top inclined plane surface) of the 2(n — l)-polyhedron is also 
smaller than 41; The shape of the result domain depends on the type to which Tj belongs. 
We shall denote the result domains by Sj. The domain Tj can always be deformed into 
the domain Sj by a transformation of the form (|2.1[) . as defined in Section 2. In the case 
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of domains of the first type, let the intersection point of n — 1 orthogonal plane surfaces be 
the pole of a system of pole coordinates r, Ox, O2, ■ • • , n -2, and let r — f(6x, O2, • • ■ , #n-2) 
be the equation of the boundary surface of T g , r — h(0x, O2, • • ■ , On-2) the equation of the 
inclined plane surface of the n-polyhedron of R™ -1 having a vertex of n — 1 orthogonal 
plane surfaces. Then the equations 

V x - X , V2-V2, ■■■ , v n -2 ~ V n -2, r - r 

J(Vl,V2, • ■ • ,Vn-2) 

represents a transformation of the domain Tj into the n-polyhedron S'j of R n_1 . For 
a domain of the second type, let cc„_i = h(xx,--- ,x n ~2) be the equation of top plane 
surface of the 2(n — l)-polyhedron and let x n -\ — f(xx, ■ ■ • , x n -2) be the equation of the 
boundary surface of T g . We then consider the transformation 

/ _ / _ t _ h(xx, • • • , x n - 2 ) 

Xi — Xx, , X n — 2 — X n — 2, X n — 1 x 



f(Xi, ■ ■ ■ ,X n -2) 

If we assume that the side length I of cube in the partition is sufficiently small, and therefore 
the rotation of the normal on the boundary surface is taken sufficiently small, then the 
transformations considered here evidently have precise the form (|2.ip . and the quantity 
denoted by e in (|2 . 1 [) is arbitrarily small. From Corollary to Theorem 10 of p. 423 of [6], we 
know that there exists a number S > depending on e and approaching zero with e, such 
that 



a k (Tj 



< S uniformly for all k, 



where 0^(1^) and ctk(Sj) are the fc-th Neumann eigenvalues of Tj and S'j, respectively. 
According to the argument as in the proof of Lemma 5.1 (i.e., (|5 .35[) ) . we see that 

4(Ej) = -L t(l n ^Ja k (Tj)), 4(E'j) = -±- t(l n y/akty)), 

where t(s) is given by and 4( E j) and 4( E 'j) (similar to <r of (|33D1 ) are the fc-th 

Steklov eigenvalue for the n-dimensional domains Ej = Tj x [0,l n ] and E'j = S'j x [0,Z n ], 
respectively. Recalling that the function t — t(s) is continuous and increasing for s > 1, we 
get that there exists a constant 8' > depending on e approaching zero with e, such that 



d( E p 



4m 



<s'. 



In other words, the corresponding fc-th Steklov eigenvalues for the n-dimensional domains 
Ej = Tj x [0, l n ] and E'j = S'j x [0, l n ] differ only by a factor which itself differs by a small 
amount from 1, uniformly for all k. Therefore, the same is true also for the corresponding 

numbers A^ e .{t) and AL, (t) of the eigenvalues less or equal to the bound r. 

s j 

The domain E'j is either a cylinder whose base is an n-polyhedron of K™ _1 having (n — 1) 
orthogonal plane surfaces with its largest side length smaller than 4/ or a cylinder whose 
base is a combination of such an n-polyhedron of K n_1 and an (n — l)-dimcnsional cube 
with sides smaller than 31; it follows from the estimates for E'j (cf. (|5.25[) — (I5.27I) 1 ) and 

Lemma 5.1 that if I is taken sufficiently small, the number A^.(t) from some t on satisfies 
the inequality 

a j e .(t) < ^rV^+CsrV- 2 
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where C\ , C2 are constants, to be chosen suitably. Thus, A E . (r) can be written as A E _ (r) = 
9(C3l n ~ 1 T n ~ 1 + Cil n ~ 2 T n ~ 2 ), where 9 denotes a number between —1 and +1 and 6*3,(74 
are constants independent of l,j and r. It follows that 

A f Ej (r) = r"- 1 [9C 3 (q - p)^ 1 + 9C 4 (q - p)l n ~ 2 \] . 

j=p+i 

As pointed out before, (q — p)l n ~ 2 < C; therefore, for sufficiently small I, (q — p)l n ~ l is 
arbitrarily small and we have the asymptotic relation 

A A F i T ) 

(5.67) lim V ' \ = 07(0, 

j=p+l 

where w{l) — > as / —> 0. For, we may choose the quantity I arbitrarily, and by taking 
a sufficiently small fixed I, make the factors of r" _1 in the previous equalities arbitrarily 
close to zero for sufficiently large r. Since 

(5.68) A%. (r) < A E . (r) for j = p + 1, • • • , q, 
we get 

(5.69) lim E » +1 5 (T) < lim E ^^ (t) =w(l) . 

T— >-\-OG T n T—t + OO T 71 

From dS35J|, (j5~50T) . (I5J31) . (|5^5)l . (|5T66]l , (I5J7) . (f5~68]) and (|5l9]l . we obtain 

(5.70) uW^-V- 1 !:^! < lim 4g <Tim-4£i 



< («„_ 1 (47r)-("- 1 )e*- 1 X;ir i |) +w(0. 

V 7 = 1 ' 



Letting / — > 0, we immediately see that X^=i -^i tends to the area |T e | of T e and 
lirm_>o TO (0 = 0. Therefore, (|5.70p gives 

(5.71) A(r) ~ ^-^ir+l.-V- 1 asr^+oo, 
or 

(5.72) A(r) ~ jf,^ 1 **- asr^+^. 

Remark 5.2. In the above argument, we first made the assumption that the boundary 
dV g of T e was smooth. However, the corresponding discussion and result remain essentially 
valid if dV g is composed of a finite number of (n — 2) dimensional smooth surfaces. 



6. Proofs of main results 

Lemma 6.1. Let gij and g'^ be two metric tensors on manifold M. such that 
(6.1) .</"'' • <■ i,j = l,— ,n 



THE WEYL-TYPE ASYMPTOTIC FORMULA FOR BIHARMONIC STEKLOV EIGENVALUES 33 



and 
(6.2) 



1 d 



1 8 



AT 



< e, i,j = !,-■■ ,n 



for all points in D, where D is a bounded domain in M.. Let 

A*i > A*2 > • • • > A*n > • • • > and //^ > /i 2 > • • • > rf n > ■ ■ ■ > 
6e positive eigenvalues of G and G' , respectively, where G and G' are given by 

f du dv 

(Gu.v) = / o-=r- tt- ds, for u and v in K, 
N ' ' J r dv dv 

(G'u.v)' = f o — — — — ds' for u and v in K! . 
Jr dv dv 

Then 

(6.3) (l + Me)-(" +1 )/ 2 (max{(l + eM),(l + Me)("+ 1 )/ 2 })~ 1 M fc <Mfe 

< (1 + Me) (n+1 ^ 2 (min{(l - eM), (1 + M £ )-( n+1 '/ 2 }) ~* fi k , 
for k = 1,2, 3,--- , 

where M and M are constants depending only on g, g' , "a^f ' "a^T an< ^ 

Proof. It follows from (|6.ip that there exists a positive constant M independent of e and 
depending only on , <?' y and D such that 

n n n 

(1 + eM)- 1 Y, 9 l 'Ut 3 < Y <0 < (! + E '/' / ' / .- 

for all points in D and all real numbers t\, ■ ■ ■ ,t n . Thus we have 

(1 + Me)-/ 2 /M < < (1 + M £ )"/ 2 ^, 

which implies (see p. 64-65 of [39]) that 

(1 + Me)- ,l/2 dR < dR' < (1 + Me) n/2 dR 

and 

(1 + Me)- {n+1 ^ 2 ds < ds' < (1 + Me) in+1 ^ 2 ds. 
(1 + Me)- (n+1)/2 [u,u] < [u,u]' < {l + Me) {n+1)/2 [u,u\. 

I d , i a 

, vi 3 — 

we immediately see that 



Thus 
(6.4) 
Putting 



'I'- % 



W\ dx i 



WW ij ) 



\g\ dx, 



{V\9\9 ij ), 



max \uJij\ < e and max|#jj| < e. 



Thus, for any u £ JC°(D) or u £ JC d (D), we have 



dxidxj 



E 



9u 
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SO that 



Ag'U ~ AgU 



d 2 u 
' dxi dxj 



du 

dx-j 



It follows that 



\A g ,u-A g u\ < e{M l \V 2 g u\+M 2 \Vgu\) ) 



where V 2 w| 2 is defined in an invariant ways as 



\Vlu\ 2 = V l V k u VzVfeU = g pl g k 



d 2 u 
dx k dxi 



1 Br 



d 2 u 



Pq dx r 



and Mi and Mi are constants depending only on g, g', ~§x~i ^§xT an< ^ ®' Thus, 



(6.5) 



J \A g ,u- A g u\ 2 dR< 2e 2 (m? J \V 2 g u\ 2 dR- 



Mi 



\V g u\ 2 dR 



Recall also (see, Section 2) that A^(D) > Af(D), where 



(6.6) 



inf 



(6.7) 



Af(D) = 



inf 



w£if<«(D), /^|v 9 H 2 rf«=i J D |V g u| 2 dR' 
and K (D) and K d (D) are as in Section 3. Let 

Xd |A 3 «| 2 di? 



(6 



(6.9) 



eS(D) = 



inf 



veK°(D), f D |V>| 2 di?.= l J D \V 2 g v\ 2 dR' 
Id \& g v\ 2 dR 



inf 



Clearly, 0?(D) > 6f(D). As in the proofs of Lemmas 2.1, 2.2, it is easy to prove that 
the existence of the minimizers to (|6 . T[) and (|6.9|) . respectively. Therefore, we have that 
Af(D) > and Qf{D) > (Suppose by contradiction that Af(D) = and Qf(D) = 0. 
Then A g u = in D for the corresponding minimizer u G K d {D) in both cases. By 
applying Holmgren's uniqueness theorem for the minimizer u G K d {D) in both cases, we 
immediately see that u e in D. This contradicts the assumption J D \V g u\ 2 dR = 1 or 
J D \ V 2 u\ 2 dR = 1 for the minimizer u 6 K d (D) in the corresponding cases). Combining 
these inequalities, we obtain 

/jA 9 u| 2 , for u€K°(D) 



D 



\A g ,u - A g u\ 2 dR < 2e 2 



M 2 



e°0D) asp). 



and 



|A ff /tt - A g u| 2 dR < 2e 2 



M 2 



m 2 \ r 



ef(D) A d (D) 
Thus we have that, for all u G K°(D) or u £ K d (D), 

(1-eM) / |A g u| 2 dR< / |A ff ,tt| 2 tfiZ < (1 + eM) I \A g u\ 2 dR 



D 



D 
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where M is a constant depending only g, g' , ^§^~i ~~§^ri ^§^~ an d D. That is, 
(6.10) (l-eM)(u,u)* < (u,u)'* < (1 + eM){u, u}*. 

By (JO]) and ([6~10)) we obtain that, for all u G K°(D) or u G K d (D), 

(1 + Me)-(" +1 )/ 2 [u,w] [«,«]' 



(max{(l + eAf),(l + Me)("+ 1 )/ 2 })((u,w)* + [u,u]) ~ «)'* + [w, «]' 

^ (l + Me)(" +1 )/ 2 [u, M ] 

" ( min{(l - eM), (1 + Me)-<"+ 1 )/ 2 }) ((it, «)* + [u, u}) ' 
which implies (|6 . 3|) . □ 

Remark 6.2. Let T and T be two bounded domains in R™ -1 , let T is similar to F (in the 
elementary sense of the term; the length of any line in T is to the corresponding length in 
r as h to 1), and let Too = Tx {a} and Too = T x {her}. It is easy to verify that 

Af(D) = h- 2 Af{D), 9({D) = ef(D), 

where D = Ti X [0, a], D = f X [0, ha], and Af(D) and Qf(D) are defined as in jB^J) and 
(|6.9|) . respectively. 



Lemma 6.3. Lei G and G' be the continuous linear transformations defined by 

e ~^>iL ds for u and v in or 



mi 



d 



(G'u,vY= [ g'^^ds for u and v m K°(D) or K d (D), 
J Te dv dv 

respectively. Let 

Mi ^ A*2 > ■ • • > f'k > • • • > and ^ > fi 2 > • • ■ > > • • • > 
6e i/ie positive eigenvalues of G and G' , respectively. If g < g' , then 
(6.11) Wk</4 /or fc= 1,2,3,.-. . 



Proof. Since g < p', we see that for any u G K°(D) or K d (D), 

(Gu,u) = Jr e g(f^) 2 rfs < Jr. g' (f^)' ds = (G^m)' 
(«,«) "( W , u )*+/ ree (|i) 2 d s -( u , w )*+/ re e'(|i) 2 ds " («,«)'' 

which implies (|6.11[) . □ 

Proof of Theorem 1.1. a) First, let (.M,;?) be a real analytic Riemannian manifold, 
and let the boundary <9f2 of f2 be real analytic. We divide the domain f2 into subdomains 
in the following manner. It is clear that the boundary dQ of the domain Q is the union of a 
finite number of closed pieces T±, ■ ■ ■ ,T p (without common inner point on the surface). Let 
U be a coordinate neighborhood which contains F,-, let Xi — Xi(Q) and a, = ^(vq) be the 
coordinates of a point Q in Tj and the interior Riemannian normal vq at Q, respectively. 
We define the subdomain Dj and surface TJ by 

Dj = {P\x{P) = x(Q) + £ n a(v Q ), Q G Tj, < £„ < a} 
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and 

T? = {P\x{P)=x{Q)+va{v Q ), QeTj}, 
where a is a positive constant. The closure of Dj is 

(6.12) Dj = {P\x(P) = x(Q) + £ n a(v Q ), Q GTj, < < a}. 

By the assumption, each f j, which is contained in a coordinate neighborhood, can be 
represented by equations 

(6.13) x t = - ,£ n -i) 

with real analytic functions ipi, i.e., it is the imagine of the closure Tj of an open domain 
Tj of M™ _1 . Hence, if a is sufficiently small, the definitions have a sense and the formula 

(6.14) x(P) = x(Q) + £ n a(v Q ), Q e f„ < £„ < a 

defines a real analytic homeomorphism of a neighborhood of the image of Dj in R" given 
by the coordinates x and a neighborhood Uj of the closed cylinder Fj in W 1 defined by 
— {£|(£i>"' >£n-i) € Tj, < < a}. Moreover, the domains D\,--- ,D P have no 
common inner points and the remainder Dq = Cl — U*j =1 Dj of O has a finite number of 
connected parts. Note that the boundary of Dq contains no part of dCl. 

Let us define the spaces K — K(d), JC and the transformation G as in Section 5. We 
shall investigate the asymptotic behavior of A{t) with regard to transformation G on space 
JC. Moreover, we define the function spaces 

K$ = {u 3 \u 3 e Hl{D 3 ) n H^Dj) n C°°(Dj), |=0on rj, and 

Au = on dDj - {Tj UTJ)}, 
#o° = {"0^0 G Hi(D ) n ii 2 (A>), = on d£> }, 
Kf = { Uj \uj e H 2 (D,), Uj = on T,, Uj = = on TJ}, 

(j = o,i,--- 

and the bilinear functionals 

(6.15) { U j,UjYj=( \A g u,\ 2 dR, (j = (),••• ,p), 

./£>., 

,„ , „. r . /" 9u 7 - <9it 7 - , . . . . 

(6.16) [Uj,Vj]j= g ~fo7~&7 > U = 1 >---,P), [uo,«o]=0, 

and 

(6.17) = + 0' = 0,---,p), 

where Uj,Vj S if? or ifj\ Closing if? and if^ with respect to the norm \uj\j = \J (uj,Uj)j, 
we get the Hilbert spaces JC® and JC d , (j = 0, • ■ • ,p). In the same manner as in Section 5 we 
can define the Hilbert spaces JC° and JC d , and then define the positive, completely continuous 
transformations G° , G d , G? and G d on /C°, JC d , /C? and /C- , respectively. Consequently, we 
can prove 

(6.18) A°(t) < A{t) < A d (r) for all t, 
and 

(6.19) A°(t) = 5>?(r), A d (r) = 5>?(r), 
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where A°(t), A d (r), A®(t) and A d (r) are the numbers of eigenvalues of the transformations 
G°, G d , G° and G d on K° , JC d , KPj and Kj which are greater than or equal to (1 + r)" 1 , 
respectively. 

Since [u ,wo]o = for all u Q G K$ or K§ and (G{ju , u ) Q = (G$u ,ii }o = [u ,u ]q, we 
immediately find that Gg = G d = 0, so that A%(t) = A$(t) = 0, (t > 0). Thus we need 
estimate A®(t) and Aj(t) for those domains Dj, where J r gds > 0. 

We can choose a finer subdivision of dQ by subdividing the domains Tj into smaller 
ones, e.g. by means of a cubical net in the coordinates £. According to p. 71 of [39], by 
performing a linear transformation $ of the coordinates we can choose a new coordinate 
system (77) such that 

g il (fj) = 6 il , (i,I = l,...,n), 

for one point fj € Tj, where Tj := <I>(Tj). Setting 0j = ipi o and a, = a, o $ , we see 
that 

(6.20) Xi(P) = <f>i(rii,- ■ ■ ,T) n -i) +r] n ai(v(rii,- ■ ■ ,7/ n _i)), 

for (771, • • • , ?7„-i) efj, < r\ n < a 

defines a real analytic homeomorphism from Ej to the image of Dj, where Ej = {77 = 
(»7i)' ' • i ?7n)|(?7i> • • • j »7n— 1) S Tj, < r) n < a} is a cylinder in K™ (This can also be realized 
by choosing a (Ricmannian) normal coordinates system at the point fj £ Tj for the manifold 
(M,g) (see, for example, p. 77 of [23]) such that a(u(r])) = (0, • • • , 0, 1) and by using the 
mapping (|6.20p .) If we denote the new subdomains of dfl by Tj as before, it is clear that 
we can always choose them and a (i.e., by letting a sufficiently small and further making 
a finer subdivision of dil, see p. 71 of [39]), so that, 

(6.21) \9*W)-g a m<e, M = V--,n, 



(6.22) 



1 d , r—rr ,^ 1 8 



j-(V\9W)\g a (v')) -^=^(M|/g 



ij = 1, 



< e, 



for any given e > 0, and all points 77' € £j. The inequalities (|6.2ip imply that 

n n n 

(6.23) (1 + M.e)- 1 < J2 9 U (v%ti < (1 + M;e) £ t? 

i— 1 — 1 i— 1 

for all points 77' G £j and all real numbers t\,--- ,t n , where Mj is a positive constant 
depending only on g d and Ej (cf. Lemma 6.1). This and formula (128) of [35] say that 

(6.24) (1 + eMj)- n / 2 \Tj\ < \T 3 \ < (1 + eMj) n / 2 \Tj\, 
where 



r < : / \Zg(j])di 11 ---dr] n - 1 , \Tj\ = df] ■ ■ ■ drj, 



are the Riemannian and Euclidean areas of Tj and Tj|, respectively, 

Next, we consider the Hilbert spaces K.j and K d . When transported to Ej , the underlying 
incomplete function spaces K j and K d are 

K] = {u\u e Hl{E ) n H 2 (Ej) n C°°(Ej), q^ = on Tf, and Am = on dEj - (Tj U Tj 7 )} 
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and 

Kf = {u \ U] e H^E,), Uj = on T h Uj = ^ = on T/}, 
respectively. The inner product, which is similar to Section 5, is defined by 

(u,v)j = (A g u){A g v)y/^jjjdr)i ■ ■ ■ dri n + g |^ |^ ^f^q) dr) X ■ ■ ■ jy„_i 

and the transformations G° and G^ are defined by 

(G^u,v)j=J g^ ^y/g(r)) dr}! ■ ■ ■ r/n-u for w,u in /C°, 



and 



(GjU,v)j= / Q——\fg(j])dr)x---na-\-, for u,t> in 



respectively. 
Put 



(6.25) g . = inf g and = sup g, 

and let us introduce the inner products 

(u,v}.= f (Au)(A g v)dr] 1 ■ ■ ■ dr] n + f g |^ dru ■ ■ ■ drj., 
J J Ei J Ti -J dv dv 



•n-l 



and 



{u,v}-= / (Au)(A g v)dr)i ■ ■ ■ dr/ n + / g — — <% ■ ■ ■ drj. 



du dv 



in the spaces K® and Kj, respectively. By closing these spaces in the corresponding 
norms, we get Hilbert spaces and K,j. Furthermore, we obtain the positive, completely 
continuous transformations G° and G^ on and Kj , which are given by 

(6.26) (G%, v) 3 ■= [ g.^^dr]i---dr] n -i, for u and v in JC ( - 
and 

(6.27) (G^u,v)j = f 7j,— ^ dm ■ ■ ■ dri n -\, for u and v in /cf, 

JTj ov dv J 

respectively. 

Let (Ufe(Gj) be the fc-th positive eigenvalue of G° and so on. According to Lemma 6.1 
and Remark 6.2, Af(Dj) and Qf(Dj) have uniformly positive lower bound when repeated 
taking finer division of D (In fact, by repeated halving the side length of every rectangular 
parallelepiped in the partition net of the coordinates r\ for each cylinder Ej, we see that 
Af(Dj) will tend to +oo, and that Qf(Dj) will have a positive lower bound). This implies 
that the corresponding positive constants Mj and Mj have uniformly upper bound when we 
further divide the domain D into finer a division, where Mj is defined as before, and Mj is 
a constant independent of e and depending only on <?, d Q Z ™ , a 9 anc j Ej as in Lemma 6.1. 
Denote by Cj (e) the maximum value of (1 + eMj) (ri+1)/2 ( max{(l +eMj), (1 +eif J ) (n+1)/2 }) 
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and (1 + eMjO (n+1)/2 (min{(l - eMjO, (1 + eM,-)" (n+1)/2 }) • Obviously, c,-(e) -> 1 as 
e — > 0. By virtue of (|6.21|) and (|6.22j) . it follows from Lemmas 6.1 and 6.3 that 

(6.28) Hk(Gj) < Cj(e) fik(Gj) 
and 

(6.29) /i fc (G$) > Ci (e)-V(G°), 
so that 

A^(r)<^(c,(e)r + c,(e)-l) 

and 

(6.30) A°(r) > ^(^(eJ-V + c^e)- 1 - l) 

where Aj (t) and 4?( T ) are the numbers of eigenvalues of the transformation (1 + r) _1 , 
respectively, and G° which are greater than or equal to By (16.181) and (|6 . 10[) , we obtain 

(6.31) £4%( e r v + c ^r l - !) < A ( r ) 

<£^%(e)T + c,-(e)-l). 

Finally, we shall apply the results of Section 5 to estimate ^(t) and A^r). Note that 

(6.32) lj(r) < Ajir) for all r > 0, 

where At is defined similarly to ([STTSTl — (|5"T7|l . It follows from ([5^1) . ([5^5]) . ([53^ and 
(j09)) - (l5TT|) that 

(6.33) !im > c n _ 1 (47r)-("- 1 )|T J |^- 1 
and 

(6.34) Hm < ^-i(4tt)-("- 1 )|T :? -|^- 1 , 

t— !-+oo T J 

where |Tj| is the Euclidean area of T 3 . By (|6~3l"]) . (|6~32"1) . flOU), ([6734]) . (|Q8]) . (|Q9l) and 
(|BT24l) . we find that 

IWAMr-f"- 1 ) < w „_ 1 (4^-("- 1 )a J ( e )^^- 1 |r J |, 

and 

hm^r-^ 1 ) >c n _ 1 (4^-("- 1 )£ J ( e )- 1 ^p»- 1 |r j |, 

T— VOO 

where Cj(e) = (1 + eMj) n ' 2 Cj(e) n ^ 1 . Note that g is Riemannian integrable since it is non- 
negative bounded measurable function on T g . Therefore, letting e — > 0, we obtain the 
desired result that 



(6.35) A(t) ~ -rf^r — / e "-Ms as t — )■ +oo. 

b) Next, since a C 2 -smooth metric 5 can be approximated in G 2 by a metric g' e which 
is C 2 -smooth on M. and piecewise real analytic (i.e., g' e is C 2 -smooth and g' e is composed 
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of a finite number of real analytic functions) in any compact submanifold of (M., g) such 
that 



gf - g a \ < e, i,l = !,••• ,n, 




< e, i, I = !,■■■ ,n, 



for all points in Z), with any given e > 0. In addition, any bounded domain D with C - 
smooth boundary can also be approximated (see, the definition in Section 2) by domain D' c 
with C 2 -smooth and piecewise real analytic boundary. Thus, the methods of Lemma 6.1 
and a) still work in this case, so that we can estimate the eigenvalues for g'* 1 in D' e . But 
for these eigenvalues (I6.35|) is true. Therefore, letting e — > and noticing that ds' t — > ds, 
we get that (|6.35[) also holds for the C 2 -smooth metric g d and D. □ 

c) With the same arguments as in the case b), we immediately see that the formula (|1.8[) 
is still true for a bounded domain with a piecewise C 2 -smooth boundary in a C^-smooth 
Riemannian manifold. 

Remark 6.4. Our method in the proof of Theorem 1.1 is new and significantly different 
from that of [3 9) . In |39j Sandgren used a technique of Lipschitz image of a convex subset 
for the harmonic Steklov problem. In our proof, Tj needn't be the imagine of a convex 
subset. Next, in order to estimate A d (r), we introduce a new counting function A? (r) as 
done in Section 5. In addition, we use the uniform boundedness of the constants Mj and 
Mj to estimate the asymptotic behavior for any finer division according to Lemma 6.1. 

Proof of Corollary 1.2. By ()1.8|) . we have 

\n - 1 

(6.36) A(X k ) ~ (vol(Sfi)) , as k +oo. 

Since A(Xk) — k, we obtain (jl.9|) . which completes the proof. □ 
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